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^ I Abstract 

Given a permutation 7r chosen uniformly from S n , we explore the joint distribution 
of 7r(l) and the number of descents in ir. We obtain a formula for the number of 
permutations with Des(7r) = d and 7r(l) = fc, and use it to show that if Des(7r) is fixed 
' at d, then the expected value of 7r(l) is d+ 1. We go on to derive generating functions 

for the joint distribution, show that it is unimodal if viewed correctly, and show that 
when d is small the distribution of 7r(l) among the permutations with d descents is 
| approximately geometric. Applications to Stein's method and the Neggers-Stanley 

■ problem are presented. 

i 

a ■ 

1 Introduction 



Consider S n to be the set of all bijections from {1, 2, . . . , n} to itself. We will often identify 
a permutation ir with the sequence 7r(l), 7r(2), . . . , 7r(n). So for instance if tt(1) = k and 
7r(n) = £, we say that ir "begins with" k and "ends with 



A permutation ir is said to have a descent at i if 7r(z) > 7r(z + 1). That is to say, if we graph 
the points (i,iv(i)) and connect them left to right, descents are the positions at which the 
connecting segments have negative slope. Let Des(7r) be the number of descents in 7T, and 
define 

(1) ("\ :=#{7re<S n :Des(7r) = d}. 

These are known as the Eulerian numbers, and have been widely studied; see, for example, 
|(;KP94I p. 267] and j( ^ar59j . 
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Bayer and Diaconis JBD92J showed that the probability that a particular permutation of a 



deck of cards occurs after any number of riffle shuffles is determined by the number of descents 
the permutation has. In |CVarj , Viswanath and the author began working to generalize that 
result to decks containing repeated cards. At one point we had occasion to consider the 
number of permutations of n letters which have d descents and begin with k. That is, 



(2) (^j :=# {vr G S n : Des(vr) = d and vr(f ) = £;}. 



The current work is an investigation of the numbers defined in Equation (J2J). We derive a 
formula in terms of binary coefficients: 



Theorem [T] // 1 < k < n, 



:)^£(-H-> io+iri 



where 0° is interpreted as 1 



which is similar to a well-known formula for the Eulerian numbers. We use the formula to 
understand how the two statistics Des(7r) and 7r(l) interact. 

If we are constructing a permutation with d descents from left to right, and d is small, a 
conservative strategy would seem to be to start with a low number, since starting with a high 
number means we will use up one of our descents near the beginning of the permutation. So 
in other words, we expect that if d is small then there are more permutations with d descents 
starting with low numbers than starting with high numbers. Similarly, if d is close to n, our 
intuition is that that starting with a high number leaves us more possibilities later on. This 
intuition turns into a surprisingly simple result: 

Theorem Q5] If it is chosen uniformly from among those permutations of n that have d 
descents, the expected value ofn(l) is d+1 and the expected value of ir(n) is n — d. 

And in Theorem [7| we find, as expected, that the sequence 




is weakly decreasing when d is small and weakly increasing when d is large. Consequently 
that sequence is an interpolation between its endpoints, which are two Eulerian numbers: 
and (j~r)- Experimental evidence (see Section [TUJ) suggests that it is a good interpo- 
lation, at least when d is close to (n — l)/2, in the sense that a normal approximation to 
the Eulerian numbers also seems to provide a good approximation to the refined Eulerian 
numbers. However, the normal approximation is good for neither set when d is small or d 
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is close to n. Theorem [7| shows that in those cases the distribution of 7r(l) is approximately 
geometric. 



The application which led directly to the current work is presented in Section Fulman 
shows in |Ful04j that certain statistics on permutations, one of which is descents, are approx- 
imately normally distributed. The main tool he uses is Stein's method, due to Charles Stein 
in |Ste86j . The thrust behind the method is to introduce a little extra randomness to a given 
random variable to get a new one. If certain symmetries are present, the result is an "ex- 
changeable pair" of random variables, meaning, essentially, that the Markov process which 
takes one to the other is reversible. Then Stein's theorems (and more recent refinements of 
them) can be applied to bound the distance between the original variable's distribution and 
the standard normal distribution. 

Fulman uses a "random to end" operation to add randomness to permutations. That is, he 
starts with a uniformly distributed permutation tt and sets 

7r' = (1,1 + 1,..., n)ir 

where I is selected uniformly from {1,2, ... ,n}. While (tt, tt') is not an exchangeable pair, 
it turns out that (Des(7r), Des(7r')) is, and this leads to a central limit theorem for descents, 
and for a whole class of statistics. 

We tried a different method of adding randomness to tt, namely, following tt by a uniformly 
selected transposition. That calculation (which is presented in Section |5J) led directly to 
Theorem El 

The Neggers-Stanley Conjecture, now proved false in general f |Bra04| IStearj ). was that the 
generating function for descents among the linear extensions of any poset has only real zeroes. 
Since a function with positive coefficients can have no positive zeroes, any combinatorial 
generating function with all real zeroes can be written in the form 

a(x + c±)(x + c 2 ) • • • (x + c n ) 

for non-negative constants a, c±, C2, ■ ■ ■ , c n . The implication, then, is that if D is the number 
of descents in a uniformly selected linear extension of a poset for which the Neggers-Stanley 
conjecture is true, then D can be written as the sum of independent Bernoulli variables. 

In Section we present several generating functions for the refined Eulerian numbers. The 
set of permutations of n which begin with k is the same as the set of linear extensions of the 
poset defined on {1, 2, . . . , n} by k < a for all a other than k. So we can find the Neggers- 
Stanley generating function for this poset explicitly, and we show that it does indeed have 
only real zeroes. We go on to show that several similar posets also satisfy the conjecture. 
(All of the posets considered were known to satisfy the conjecture by theorems of Simion 
|Sim84| and Wagner |Wag92| .) 
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2 Basic Properties 



If 7r(l) = 1, then 7r(l) is certainly less than 7r(2), so all descents are among the final n — 1 
numbers. And if 7r(l) = n, there is certain to be a descent between 7r(l) and 7r(2). So we 
know some boundary values: 

. . /n\ /n — l\ , /n\ ln — \ 

for n > 1. Also, it is immediate that 

w EC) =(»-i)' 



(5) £ 



7i\ / n 
d/u = \d 



Let p & S n he the reversal permutation: p(i) — n + l — i. Then p7r is the same as n but with 
i replaced by n + 1 — i everywhere. As a result, pir has a descent wherever 7r has an ascent, 
and an ascent wherever 7r has a descent. So Des(/wr) = n — 1 — Des(7r). Since 7r i— > p7r is a 
bijection from S n to itself, it follows that 

<«) (:>-(.-"_- 

Note we could have obtained the same result from the map n \— > 7rp, since reversing n changes 
ascents to descents and also reflects their positions about the center. 

Let 

(7) (^j :=#{neS n : Des(7r) = d and n(n) = k} . 

Both transformations yield symmetric identities for the refined Eulerian numbers. If 



then 

p7r(l) = n + 1 — k and Des(p7r) = n — 1 — d 

n — 1 — d 



7T(1) = 


= k 


and 


n + 1 


- k 


and 


Tip{n) 


= k 


and 


n + l 


-k 


and 


n 






■1-d, 




-k 



from which it follows that 

/n\ I n \ I n \ /n\ 

(8) 



d/ k \n-l-d/ n+1 _ k \n-l-d) \d) 
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3 Recurrences 



Assume n > 1. Let 

Tk '■= {vr G £„ : tt(1) = k and Des(7r) = d} 
Tk,e '■= {vr G iS n : tt(1) = fc, 7r(2) = £, and Des(7r) = d} 

and let 7r G Tj.^. If £ < k, then there is a descent between 7r(l) and vr(2), so there must be 
d — 1 descents in the "tail", 7r(2),7r(3), . . . , 7r(n). The tail begins with £, which is the £th 
largest value in the tail, so we must have 

#Tm = (d - 1 

when £ < k. Likewise, if £ > k, there is no descent between 7r(l) and vr(2), so there must be 
d descents in the tail. This time £ is the {I — l)st largest value in the tail, so 

#T M = ( n 1 

" / £-1 

when £ > k. Of course T fc is the disjoint union of the T fc /, so 



or, more succinctly, 

(») C)=E 



\ n— 1 , , 

n\ \ ^ / n — 1 



d/k f^\d-[£<k] 

where the bracket notation follows |GKP94j : [A] is 1 if A is true and if A is false. (Knuth 
refers to this as Iverson notation in |Knu92j . and traces its origin to |Ive62j .) Note Equa- 
tion Q fails when k < 1 or k > n, in which case (^/ k = 0. 

Now suppose 1 < k < n — 1 and 7r £ 5 n begins with fc. Swapping k with in the sequence 
7r(l), 7r(2), . . . , 7r(n) preserves descents for most 7r; the only exception is when it (2) — k + 1, 
in which case a new descent is created. If we eliminate that case, the swap map is a bijection 
from T k \ T ktk+ i to T fe+1 \ T fe+ljfc , as those sets are defined above. Substituting sizes for sets, 
we have 

/n\ /n — l\ / n \ /n—V 
(10) W*"\ d ) k = \d/ k+1 ~ \d-l /k 

Equation (jl(jp is valid as long as k ^ and k ^ n. (If k < or k > n, all terms are 0.) 

A well-know recurrence for comes from considering what happens when you insert n into 
an element of if> n -i : 
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We can get a similar recurrence for the refined Eulerian numbers by considering what happens 
when you insert n into an element of S n -i which begins with k: 

/n\ , , — 1\ . /n — 1' 

In other words, one way to get an element of S n which begins with k and has d descents is 
to take an element of SVi-i which begins with k and has d descents, and insert n at a descent 
or at the end (d+1 choices). The other way is to start with an element of S n -x which begins 
with k and has d — 1 descents, and insert n at an ascent (n — d — 1 choices). Equation (|T2*j) 
fails when k = n, since a permutation of S n -i cannot begin with n. It is valid for all other 
values of k. 



4 Formulas and Moments 



There is an explicit formula for the Eulerian numbers in terms of binomial coefficients: 

See for example, |&KP94 ( p. 269]. [Aside: Equation (J13)) follows from Equation (Jll)) . which 
means that it is valid for all values of <i, even if d < or d > n]. So we have 



(15) /rA /n ~ 1 



n-1 



d \d 



These suggest a formula for (^) fc : 



Theorem 1 // 1 < k < n, 



tie) C^D-T^^n'+ir 



where zs interpreted as 1. 

Proof. Fix fc > 1. The theorem is true for n = k by Equation (|15jl. Suppose it is true for 
some n > k. Then 



- {,> - (,] [ d - n i-j) +{d+1 \ d -3 
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The quantity in brackets reduces to (j + 1)Q^), so the theorem is true for all n > k by 
induction. □ 

Note that we assumed nothing about d; Equation (fTBj) is valid even if d < or d > n. 

From Equation (fTBj) we can deduce a formula for the mth "rising moment" of 7r(l) when 
Des(7r) is fixed. Assume 7r is chosen uniformly from S n , and let 

(17) fx m := E Dcs(7r)=d 7r(ir 

where x m = x(x + i)(x + 2) • - • (x + m — 1). 

Lemma 2 



:«) (>-^E(-^( d !,)|( m r)^ 



Proof. From Equation 



(the last by setting r = k - 1). But (j + I)™- 1 -'' = { n ~l' r )j s . So let £ = r + s and 

we have 



Let be a north/east lattice path from (0,0) to (m + n — £) (see Figure Q). The 
number of such paths is ( m ^ n ) • If r is the height at which <p crosses the line x — m + |, then 
cf) consists of a path from (0, 0) to (m, r), a horizontal segment, and a path from (m + 1, r) 
to (m + n — £, £). Counting the possibilities for the parts yields the identity 

Substituting Equation (J20|) into Equation (fT9| yields the desired result. □ 
Note that the last sum in Equation (|T8|) is a truncated binomial expansion of (j + l) m+n . 

Theorem 3 If it is chosen uniformly from among those permutations of n that have d 
descents, the expected value ofn(l) is d + 1 and the expected value of it(n) is n — d. 
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m 



x = m + 2 



J 



x 



m + 1 m + n — 



Figure 1: A north-east lattice path from (0, 0) to (m + n — £, £). (All edges are either north 
or east.) 



Proof. The expected value of 7r(l) is /^i, and 

n 

s.d-3 



i>o ^ J ' e=o 



i>0 \ J/ j> 

The term for % = is 0, so let j = i — 1 and combine 



= X>1)^(J + 1)" [ ( d ! •) (i + 1) + ( rf _ " _ J (n + j + 2) 

The quantity in brackets simplifies to (d+ 1)(™^), so 

n + 1 
d-j 

/Xi = E Des W =d 7T(l) = d + 1. 



(J ft =(j+i)D-i)^ + ir(" + !) = (-+D'" 



i>o 



Therefore 



For the second part, 



= E DcsW=n - 1 - d 7r(l) = n - d. 



□ 



5 Application Using Stein's Method 

Charles Stein developed a method for showing that the distribution of a random variable W 
which meets certain criteria is approximately standard normal. His technique has come to 
be known as Stein's method; sec [Stc86] or |DH04j for more explanation than can be given 
here. 

In its most straightforward form, Stein's method requires finding a "companion" random 
variable W* such that (W, W*) is an exchangeable pair, meaning that 

(21) F(W = w,W* = w*) =F(W = w\W* =w) 

for all values of w and w*. If we can find such a W* and if, in addition, there is a A between 
and 1 such that 

(22) E W W* = (1 - X)W 

(that is, the expected value of W* when W is fixed at some value is 1 — A times that value), 
then we may apply Stein's method. 

We are interested in showing that if tx is chosen uniformly from S n , then the random variable 
D = Des(7r) is approximately normal. This has been proven before, and in more generality; 
see |Ful04j for references. We will demonstrate the set-up for Stein's method — that is, finding 
a companion variable and showing that it satisfies Equation (j2"Tj) and Equation (|2*2*|) . From 
there, applying the method would proceed as in |Ful04| . 

Often the companion variable in Stein's method is defined by adding a little bit of randomness 
to the variable we are interested in. In this case, let r be selected uniformly from among the 
transpositions in S n , independently of tx. Then tix is uniformly distributed over S n , and for 
any u,v E S n , 

P(tt = u, tix = v) = P(tt = u,t = vu~ l ) = P(vr = u)P(r = vu~ l ) 
P(tt = v, tix = it) = P(tt = v,t = uv' 1 ) = P(tt = u)P(r = uv _1 ). 

Both right-hand sides are (ri!) _1 (™) 1 if vu~ x is a transposition and otherwise, so (ix, tix) 
is an exchangeable pair. 

Let D* := Des(r7r). Since (7r,T7r) is an exchangeable pair, (F(tx), F(ttx)) is exchangeable for 
any function F. So (D, D*) is exchangeable. For 1 < i < n — 1 let 

Di = [ir(i) > 7x(i + 1)] and D* = [nr(i) > ttx(i + 1)] 

be Bernoulli random variables; then D = Y17=i and D* = Y^i=i D*- 
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Fix it and i and let a = 7r(i), b = ir(i + 1). If a < b, the only ways for rir(i) to be bigger 
than TTr(i + 1) are if r swaps a with something bigger than b (n — b ways), if r swaps b with 
something smaller than a (a — 1 ways), or if r swaps a with b. So 

E D --"(£>* - A) = V(D- = 1|A = 0) = ILt'W " + 



G) 



and similarly if a > 6 

E A=1 (A - A) = -P(A* = 0|A = i 

So in general 



A=lm* inw-v „m i> U + 7r(i + 1) - 7r(i) 







,D ifT * ^ 7T(i) - 7T(i + 1) , 2(1 -2A 



E^(A - A) = w / n \ 7 + 



G) ' 

Summing now over i causes the terms to telescope: 

W { D> - D) = £ W(D> - A) = ^=f^ +2 - 
which allows us to apply Theorem El 

YPgr - d) = k°f(p- - o) = E °" (1) ~ E ^ (n) + 2 - 4D 



« - 1 



n(n — 1) n — 1 n 

The mean and variance of Des(7r) are \l :— (n — l)/2 and a 2 := (n + 1)/12 respectively, so 
the variables 

W := Des(7r) " M and := 

a cr 

have mean and variance 1. Then (W, W 7 *) is an exchangeable pair and 

E W=W {W* -W)= E D=aw+ " ( D *~V _ ) = Ie^=™+/^(£)* _ £)) 

\ a a J a 



which is to say 



an n 



So if W* is obtained using the "random transposition" method described here, (W, W*) will 
be an exchangeable pair satisfying Equation (|22j) with A = 4/n. One can now proceed with 
Stein's method and show that W is close to being a standard normal random variable. 
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6 Generating Functions 



It follows from Equation (JT5)l that 



?? 



x 



d+l 



1 •./•r : ^. / ".r' 



and therefore that 

: = J2 a n(x)z n /n\ = ^(l-i) n+1 ^jV2 n /n! 

n>l n>l j>0 

= (1 - x) ^ - x)z) n /n! = (1 - x) J^ xj {e 3(1 - x)z - l) 

j>0 n>l j>0 



(1-x) 



1 - xe^- x > 1 



1 — x 



1 - xe^ 1 "^ 2 



x — e~^-~ x > z 

There is some disagreement in the literature about what a should be. We have avoided the 
problem by not including it in the sum. 

There are various ways to define generating functions for the (^} k , depending on which 
variables are kept constant. 



Theorem 4 

(23) 



e - \ i 

n,d,k N ' k J0 

where 9 = exp j ^-7^ J £ j • 



tit 



t i-i/» 



Proof. Let -B(x, y, z) be the left-hand side of Equation (J23J) . Note the sum is over all integers 
n, d, and k. So 



(24) 



(iT 1 - 1)^ + (i - *)B 



E 

n,d,k L 



71 + 1 

d 



n + 1 



+ i \ d / k \d/ k \d 1/ ,_ 



x d y k z n /n\ 



Let 5"(n, cf, fc) be the bracketed quantity. It is clearly if n < 0, and if n = 0, 

1 if d = 0, k = 
S^djJfe) = <{ -1 if d = 0,Jfe=l 
otherwise 
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so n = contributes 1 — y to the sum on the right-hand side of Equation (jiMj) . If n > 1, 
then by Equation (jlOj) . S^n, <i, /c) is unless fc = Oorfc = n+ l,in which case 



S(n, d, 0) 
Therefore 



71 + 1 

d 



\ and S(n, d, n + 1) = — / 
a/ \ a 



95 



(r 1 -l)^ + (l-^)5 = l-Z/+ E 



n 



x z 



n>l,d n>l,d 

1 — y + x _1 A(x, z) — yA(x, yz) 



n+l 



n 



n 

d- 1 



d-1 



x d y n+l z n /n\ 



l + x 



[1-x) 



x — e~^~ x ^ z 



y 



i + 



l 



.re 



-(l-x)yz 



X 



x — e~^~ x ^ vz 



x - e-^-^y* x - e -^- x > 



Let a = -=T3f - Then 0, as defined in the theorem, is e a2 . Dividing by y x — 1 and multiplying 
through by 9 gives 



which is to say that 



dB 

9— + a9B = a9 

oz 



ye 



-(l-x)yz 



1 



2/0 







x - 9y~ 1 x - 9 1 - 1 /y 

Differentiating the integral on the right-hand side of Equation 



9" 



dt 



d&y 



d_ 

dz J a x — t x ~ x ly dz 

ay9 





89 


1 


Vv 


dz 


x - 9 l ~ l /y _ 



a9 



d 



x - 9y~ x x - 9 1 - 1 /y dz 

Since 9B and the integral have the same derivative with respect to z, and they both vanish 
when z = 0, they are equal. □ 

Here are three more generating functions. They can all be found by plugging in Equation 
and switching summation signs. 



(25) 
(26) 
(27) 



E 

d 

E 

k 

E 

d,k 



X 



k x j 



j>0 



y =y 

xv = (i-xr y 'E u + 1)n r Uy)r x> 



j>0 



J + i - jy 



12 



We can now prove a special case of the Neggers-Stanley conjecture. Define the descent 
polynomial of A C S n to be 

F A (x) = ^x DesW . 

7TgA 

Let P be a poset of n elements with labels 1,2, ... ,n. A linear extension of P is an ordering 
of 1, 2, . . . , n which preserves the ordering of P; that is, a tc 6 5* n which is such that if z <p j 
then i appears before j in the list 7r(l), 7r(2), . . . , 7r(n). If £(P) denotes the set of linear 
extentions of P, then Neggers and Stanley |StaOO| p. 311] conjectured that for any poset, 
every zero of Fc{p) is real. 

The conjecture has been shown to be false in general [Bra04. IStearj . But we can prove it is 
true in a certain special case. 

Theorem 5 If P n ^k is the poset with Hasse diagram 

1 2 ••• k - 1 k + 1 k + 2 ■■■ n 




k 



then Fc(p nk ) has only distinct real roots. 



Proof. For u, v > let 

d 



7TtO U -)-,;-)-l 

7r(l)=«+l 



Then setting u = k — 1, v = n — k yields the polynomial in question. If v = 0, c U)V counts 
the reversal permutation p, which has (u + v + 1) — 1 = u descents. Otherwise, if v > 0, c UjV 
doesn't count p but it does count the permutation 

u+ l,u, u — 1, . . . , 1, u + v + l,u + v, . . . ,u + 2 

which has u+v-1 descents. So 



deg(c n) „) 



u if v = 

u + w — 1 iff>0. 



Similarly, if it = 0, c U)V counts the identity permutation, which has no descents. Otherwise 
it doesn't count the identity but it does count 

u + l,u + 2, ... ,u + v + 1,1,2, ... ,u 

which has 1 descent. So x \ c 0)V (x) and if u > 0, x \ c UjV (x) but x 2 \ c UjV (x). Now let 
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Dx ( 



h ,o( x ) 


= (1 


— x) 


-1 




h 0> i(x) 


= (1 


— x) 


-2 




h 0>2 (x) 


= (1 


— x) 


- 3 ( 


[1 + x) 




= (1 


— x) 


-4, 


1 + Ax + x 2 ) 




= (1 


— x) 


" 5 ( 


[8x + Ux 2 + 2x 3 ) 




= (1 


— x) 


- 6 ( 


[8x + QOx 2 + A8x 3 + Ax 4 ) 


h3,3(x) 


= (1 


— x) 


- 7 ( 


'8x + 160x 2 + 384x 3 + 160x 4 + 8x 5 ) 



Figure 2: The construction of ^3,3 (x) as described in the proof of Theorem 0D The zeroes of 
each function are plotted on the right, using an inverse tangent scale. Since each function 
is generated from the previous one by applying either the Dx or the xD operator, Rolle's 
Theorem guarantees that the zeroes must interleave. By a counting argument, all the zeroes 
of each function must be real. 

Note that c Uj „(l) = # {it G S u+v+ i : 7r(l) = u + 1} = {u + 1>)!, so c U)V does not have a zero 
at x = 1. Therefore h UtV has exactly the same zeroes as c UtV , plus a pole at a; = 1. By 
Equation (J25J), 

K, v {x) = Y,3 u {3 + l) v x 3 - 

If D represents differentiation with respect to x, we have 

(xD)h UjV (x) = h u+hv (x) and {Dx)h u>v (x) = h UtV+1 (x) 

and so 

ho tV {x) = (Dx) v hofl(x) and h u<v (x) = (xD) u h 0tV (x). 

hofl( x ) — (1 — x)" 1 an d h 0t i(x) = (1 — x)~ 2 both have no zeroes. Suppose v > 1 and h 0v 
has only distinct real zeroes. Since deg(c 0il) ) = v — 1 and x \ Cq iV (x), xc 0)V (x) and xh 0tV (x) 
have v distinct real zeroes. By Rolle's Theorem, (Dx)ho tV must have v — 1 distinct zeroes 
interlaced between those of ). Furthermore, the denominator of xh ^(x) has degree 

v + 1, so x/io,t;(x) approaches as x — ► 00. Therefore its graph must turn back toward the 
x-axis somewhere to the left of its leftmost zero, at which place there must be another zero 
of (Dx)ho tV . So we have found v real zeroes of ho jV+ i, and that accounts for all its zeroes. 

Applying the xD operator goes similarly. Given that h UjV has d distinct real zeroes, by 
Rolle's Theorem Dh UjV (x) has d—1 interlaced zeroes. Since the numerator of h UtV has degree 
smaller than the denominator, h u>v must turn back toward the axis to the left of its leftmost 
zero, which accounts for one more zero of Dh u v . Finally, (xD)h UyV has one more zero at 
(which is distinct from the others since x 2 \ h u>v and therefore x \ Dh UjV ). So we have found 
d + 1 real zeroes of h u+ i tV , and that accounts for all of the zeroes. □ 



Corollary 6 The same can be said for the poset 
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Proof. The result of turning a poset upside-down is to reverse all its linear extensions, 
which changes ascents to descents and vice-versa. So if F(x) is the descent polynomial of 
the original poset, the descent polynomial of the new poset is x n ~ 1 F(x~ 1 ). So the roots of 
the new polynomial are the inverses of the roots of the original. □ 



7 General Behavior 

We can say in general how the sequence (^) , (^) • • • , ( r ^) 1 behaves. 

The set of numbers (^), , for n fixed, is very nearly unimodal if arranged appropriately. 



Theorem 7 Fix n and d. Then 



(i) 


lfd = 0, 




o = < 


/n\ 
\d/n 




• = c> 2 < a = i 


(ii) 


Ifl<d<(n- 3)/2, 




(d)n 


<( 


d/n-1 


<---<a 


(iii) 


If n is even and d = (n - 


-2)/2, 




< ■ 


••<( 


d)2 ~ (d)l 


(iv) 


If n is odd and d = (n — 


i)A 


/n\" 
\d/n 


< • 


••<( 


d)(n+l)/2 > " " > \d)i 


(v) 


If n is even and d = n/2 




(d)n 


= ( 


n\ 

dl 7i—l 


>•••><»! 


(vi) 


J/(n + l)/2 <d<n-2 




(d)n 


>( 


n\ 

dl n—1 


>-->a 


(vii) 


If d = n — 1, 




1 = < 


/n\ 
Wn 


>(2 


u = - = a = o- 



Proof, (i) follows from the fact that the identity is the only permutation with descents, 
(v), (vi), and (vii) follow from (iii), (ii), and (i) respectively because (^) fc = ( n _"_ d ) n x _ k - 

Let f n (x) = , which means that f n (nd — k) = \£) k ii < d < n — 1 and 

1 < k < n. Figure El shows the graphs of fe(x) and f^ix). Each monochromatic section is 
a sequence of the form (^) , (^) _ 1( • • • , ( T d) 1 - Note the graphs plateau where one sequence 
meets the next. Since (^V = ("l^ 1 ) = (^x) , each sequence begins where the previous one 
ends. The content of the theorem is that /„ is basically unimodal. That is, the sequences 
on the left increase, those on the right decrease, and those in the middle behave according 
to (iii) through (v). 

The theorem is true for small n by inspection. By Equation JSJ), 
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66 




d=l d=2 



384 



il 



y =/,(*) 



-d=2,k=5 



Figure 3: The graphs of /e(^) and fi{x), where f n (nd— k) = (^),, as defined in Theorem[7| 



Let i = £ + n — km the first sum and £ — k in the second and we have 

re— 1 re— fc re— 1 



n + 1 
d 



E /M*- 1 ) 



(z — n + 



*)) + E a™* - (* + ^ = E ^ - fc 



i=re— fc+1 



8=0 



i=0 



So imagine a caterpillar of length n crawling on the graph of y — f n {%), as shown in the top 
graph of Figure El If his head is at x-position nd — k, the equation above says that the sum 
of the heights of his segments (or his total potential energy) is ("^ ) . If he were to take 
a step forward, his total energy would be That would be an increase in energy if 

the new height of his head is higher than the current height of his tail. The theorem now 
follows easily by induction. □ 



8 Behavior if d <C n 



If d is much less than n, and it is selected at random from those permutations of n letters 
which have d descents, then the distribution of 7r(l) approaches a geometric distribution 
uniformly, in the following sense. 



Theorem 8 Fix an integer d > 0. Suppose ir n is chosen uniformly from those permutations 
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of n letters which have d descents. Then for any e > there is an N such that 

Wl) = k) 



28 



1 



< e 



(1 — p)p k 1 

for all integers n and k with n> N and 1 < k < n, where p = 



Proof. For < j < d, let Pj(n) = {-l) d ~ j . Then by Equation JEJ and Equation $F$, 

I r>.\ , , , < — - 

(29) 
(30) 



d k -\d+\) n ~ k p j 

0<j<d 



n 



0<j<d 



" ]{ 2) [l+i 



Since (1 - p)p k 1 = d k 1 /(d + l) k , the left-hand side of Equation (j2BJl is 



, n— fc 



£o<;< d ^ + i)(fii) n 

and since Pdiji) = (™) = 1, the last term of both sums is 1. Therefore we have 

E„< 3 <, Hn) (if- 1 (i±l)-* - />,(„ + 1) (i±i)"l 



i + E„ s <,^(« + i)(S)" 

Since j/d < (j + l)/(d + 1) when < j < d, that's bounded above by 



Ho<j<d 


\p 3 (n)\mr 1 + 











Now each term in each sum is a polynomial in n times a decaying exponential in n. So both 
sums go to as n goes to infinity. □ 



Corollary 9 The total variation distance between the distribution of 7r n (l) and the geometric 

d 
d+l 



distribution with parameter p = -^-r approaches as n approaches infinity 



9 If Both Ends Are Fixed 



We might now ask about the number of permutations with d descents whose first and last 
positions are fixed. Let 



# {ir G S n : Des(7r) = d, 7r(l) = k, and n(n) 
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Theorem 10 Suppose 1 < k < k + m <n. Then 



d/ k \ d I \d/ k+m \d-l 



Proof. Let if) e S n be the n-cycle (n, n — 1, . . . , 2, 1). Then for any 7r e S n , 

WJ \n if 7r(i) = 1. 

(Imagine a device like a car odometer, with a window and n wheels, on each of which are 
painted the numbers 1 through n. it can be represented by turning the ith wheel until 7r(i) 
shows through the window, for all i. If one then rolls all the wheels backward a notch, ifj7i 
shows through the window. For this reason we will refer to the transformation 7r i— > ipn as a 
rollback.) 

If 1 < i < n, let Di(ir) = [ir(i) > ir(i + 1)]. The pair n(i + 1) has one of four types: 



Type 


AW 


A(^) 


AftH- AW 


A 


1 < 


7T 


'i) < 7r(i + 1) 











B 


1 < 


7T 


+ < 7r(i) 


1 


l 





C 


1 = 


7T 


;i) < ?r(i + 1) 





l 


1 


D 


1 = 


7T 


+ < 7r(i) 


1 





-1 



Most pairs are of type A or B. 7r will have one pair of type C unless n{n) = 1 and one pair 
of type D unless 7r(l) = 1. Therefore 



n-l 



Des(^vr) - Des(7r) = ^ A(^) ~ A(vr) 



1 if tt(1) = 1 
-1 if ir(n) = 1 
otherwise. 



Let 



Pa := { n e 'S'n : tt(1) = a and 7r(n) = b} 
Q b := {tt G 5 n _! : 7r(n) = 6} . 



Consider the following sequence of bijections: 

pfc+m rollback p fc-i+ m rollback rollback p 



1+m 



rollback 



pm 

n 



shorten 



where "shortening" a permutation means removing n. (See Figure /reffig: rollback for an 
example.) The first k — 1 rollbacks all preserve Des, and the final one increments Des. But 
the shortening decrements it again, since it removes n from the front of the permutation. 
Therefore the net effect, across the whole sequence, is to preserve Des. So (^)^ +m = ("rf 1 )™ 
for all d. 
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Descents 



Descents 



3 7 9 


1 8 


4 5 


2 6 



2 6 8 9 


7 


3 4 15 


15 7 8 


6 


2 3 9 


4 




9 


4 6 7 


5 


1 2 8 


3 








4 6 7 


5 


1 2 8 


3 



8 9 13 


2 6 


5 7 


4 




7 8 9 


2 


1 5 


4 6 


3 


6 7 8 


1 9 


4 


3 5 


2 




5 6 7 9 


8 


3 


2 4 


1 



4 


5 


6 8 


7 


2 


1 3 


9 
















4 


5 


6 8 


7 


2 


1 3 





Figure 4: Examples of the actions of the bijections described in Theorem for n — 9. 
Vertical lines show the positions of descents. If 7r(l) < 7r(n), as at the top left, then the 
permutation is "rolled back" until n appears at the front, and then n is removed. In each 
of the rollbacks but the last, one of the internal bars moves one position to the right, to 
accomodate a 1 changing to an n, but the total number of descents stays the same. Only 
when the number in the first position changes from 1 to n do we gain a descent, but it vanishes 
again when we remove n in the last step. The procedure is is similar when 7r(l) > ir(n), as 
on the right, but the last rollback eliminates a descent, and removing n leaves the number 
of descents unchanged. 



The second part of the theorem follows from the bijective sequence 

pk rollback p k-i rollback rollback pl rollback p „ shorten „ 



where Q a = {tt e S n : 7i"(l) = a}. Here the final rollback decrements Des, and the shortening 

k+rn \d-l/ m ' 



leaves it unchanged. So ( 1 }) h k+m = (^-i) • ^ 



Corollary 11 // 1 < k, £ < n and k is the poset on {1, 2, . . . , n} defined by k <p a <p £ 
for all a other than k and £, then the descent polynomial of C(P^ k ) has only distinct real 
zeroes. 



Proof. If I — k + m, then the polynomial in question is 

which was shown to have real distinct zeroes in Corollary El As in that corollary, it follows 
immediately that turning the poset upside-down inverts the roots of the polynomial, leaving 
them real. □ 
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10 Remarks 



In Section El we noted that if 7r is uniformly distributed over S n then the distribution of 
D = Des(7r) is approximately normal. Thus the normal density function 



Experimental evidence for n < 200 suggests that this is in fact the case. So while the 
distribution of 7r(l) given Des(7r) is by no means normal, it does seem to behave like a 



More generally, there may be some underlying curve which the Eulerian numbers, properly 
normalized, can be said to approach as n grows large. It will look like a bell curve, but 
not be exactly normal, since the normal approximation is not very good when d -C n. If 
so, it seems likely that the refined Eulerian numbers presented in this paper can be said to 
approach points on the same curve. 
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